NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 

aRe  Used  for  any  purpose  other  than  in  connection  with  a definitely  related 

GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


Reproduced  by 

DOCUMENT  SERVICE  CENTER 

LDIN6.  DAYTON.  2.  OHIO 


^ Derivation  of  Atomic  Co-ordinates 
cri  c/j  ^ Helical  Structures 


•'■-u 


A Single-Crystal  Adaptor  for  the  Norelco 
High-Angle  Diffractometer 


Technical  Reports  71  and  11 
Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 


September,  1953 


AD  No./2-^ 

ASTIA  FILE  COPY 


Technical  Report  7 1 

Derivation  of  Atomic  Co-ordinates  in  Helical  Structures 

by 

H.  J.  Grenville-Wells 

Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 
Cambridge.  Massachusetts 


i 

» 

f 

i 


I 

r 

k 


O.N.R.  Contracts 


NSori-07801 

N5ori-07858 


f 

i 

September,  1953 


DERIVATION  OF  ATOMIC  CO-ORDINATES  IN 
HELICAL  STRUCTURES 


by 

H.  J.  Grenville-Wclls 

Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Abstract;  A general  method  is  developed  for  obtaining  the  co-ordinates  of  atoms 
in  helical  structures.  The  helix  is  treated  as  the  repetition  of  a struc- 
tural unit  by  an  n-fold  screw  axis,  where  n is  continuously  variable  in- 
stead of  being  restricted  to  integral  values.  Expressions  for  evaluating 
the  co-ordinates  in  terms  of  an  orthogonal  axial  system  are  derived,  and 
graphs  which  facilitate  the  use  of  these  expressions  are  given. 

All  crystal  structures  belonging  to  space  groups  which  contain  screw  axes 
can  be  formally  considered  as  helical  structures.  In  the  case  of  n-fold  screw 
axes  where  n is  an  integer  with  values  Z,  3,  4 or  6,  the  co-ordinates  of  the  atoms 
which  are  repeated  by  the  screw  can  be  obtained  quite  simply  in  terms  of  the 
co-ordinates  of  general  and  special  positions  for  the  appropriate  space  group. 
These  co-ordinates  are  listed  in  Vol.l  of  the  International  Tables  for  X-Ray 
Crystallography  (1952). 

In  more  complex  structures  such  as  proteins,  however,  in  which  it  is 
possible  that  there  may  be  helices  approximating  the  repetition  of  structural 
units  by  n-fold  screw  axes  where  n is  nonintegral,  the  co-ordinates  of  atoms  in 
successive  units  cannot  be  readily  listed,  since  n is  continuously  variable.  The 
co-ordinates  are  needed,  however,  for  calculations  of  bond  lengths  and  bond 
angles,  and,  since  they  vary  smoothly  with  changes  in  n,  graphs  can  be  drawn 
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which  greatly  reduce  the  amount  of  work  involved  in  obtaining  the  co-ordinates 
for  a given  structural  unit  repeated  by  an  n-fold  screw  axis  for  a particular 
choice  of  n. 

The  need  for  a systematic  derivation  of  this  type  was  recognised  by 
Dr.  Barbara  W.  Low  of  Harvard  Medical  School,  in  the  course  of  her  work  on 
proteins,  and  the  specific  form  in  which  the  derivation  is  given  below  was  de- 
veloped to  analyse  polypeptide  chain  helices. This  concept  can  be  adapted 
easily  to  deal  with  structures  containing  infinite  chains,  such  as  Te  and  Se,  or 
with  the  small  deviations  from  true  screw  symmetry  involved  in  spiral  distortion 
in  crystals  which  contain  screw  axes, such  as  quartz.  The  mathematical  deri- 
vation of  it  is  therefore  given  below. 

Consider  first  the  structural  unit  which  is  to  be  repeated  by  the  screw 
axis  (in  this  case  the  residue  * • C|  ’ > which  is  shown  inscribed 

in  a circle  of  diameter  in  Fig.  1).  The  residue  is  assumed  to  be  planar,  and 

its  plane  is  always  taken  to  be  parallel  to  the  helical  axis.  Relaxation  of  this  con- 
dition may  be  introduced  later  on  in  the  calculation. 

If  one  end,  C^,  is  kept  fixed,  the  residue  may  be  tilted  either  up  or  down. 

If  P is  the  pitch  of  the  final  coil,  i.e. , the  vertical  translation  per  complete 
turn,  and  n is  the  number  of  residues  per  complete  turn,  then,  in  a given  helix, 
this  vertical  translation  per  residue  is  P/n.  If  the  residue  is  tilted  upwards  so 
that  (solid  line)  is  its  horizontal  projection  (normal  to  the  helical  axis), 

then  the  CO  and  NH  groups  are  almost  vertical  and  parallel  to  the  helical  axis. 

A value  of  P/n  = L14  A was  used  to  construct  Fig.  1. 

When  the  residue  is  tilted  downwards,  so  that  (broken  line)  is  hori- 

zontal, then  the  CO  nnd  NH  groups  are  both  sharply  inclined  to  the  helical  axis, 

II 

1)  B.  W.  Low  and  H.  J.  Grenville-Wells,  Generalized  Mathematical  Relation- 
ships for  Polypeptide  Chain  Helices.  The  co-ordinates  of  the  ir  helix,  " 

Proc.  Nat.  Acad.  Sci.  Wash.,  1953,  in  press. 
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whereas  C^C|  and  are  nearly  parallel  to  ^^^2' 

It  should  be  noted  that  these  two  directions  of  tilt  define  two  different 
helices,  which  nevertheless  have  the  same  value  of  P/n.  Systematic  variation 
of  P/n  will  therefore  result  in  two  distinct  series  of  helices,  which  cannot  be 
transformed  into  one  another  by  winding  more  tightly  or  more  loosely. 

This  will  always  be  true  unless  the  line  is  a mirror  line  in  the  struc- 

tural unit,  in  which  case  the  upward  and  downward  tilts  associated  with  a par- 
ticular value  of  p/n  produce  the  enantiomorphous  pair  of  helices  corresponding 
to  a left-handed  and  right-handed  screw.  This  is  not  the  case  for  the  residues 
of  polypeptide  chains,  so  that  two  series,  which  we  have  called  the  "a"  and  "y" 
series  because  they  contain  Pauling's  "o"  and  "y"  helices,  respectively,  must 
be  considered. 

Atomic  co-ordinates  for  all  atoms  in  the  helix  are  derived  in  two  stages: 

(1)  the  co-ordinates  of  atoms  in  the  basic  structural  unit  are  obtained  with  the 
fixed  end  of  the  unit  taken  as  origin  (C^  in  Fig.  1);  (Z)  the  co-ordinates  of  atoms 
in  successive  units  are  obtained  by  the  operation  of  the  n-fold  screw  on  these 
"basic"  co-ordinates. 

First  stage.  If  the  basic  unit  is  placed  in  the  plane  y = 0,  only  the  x and  z 
co-ordinates  remain  to  be  determined.  (If  a nonplanar  unit  is  to  be  considered, 
this  would  no  longer  be  true.) 

It  can  be  seen  in  Fig.  1 that  for  a given  structural  unit,  the  angles  6^,  . . . 

and  the  distances  of  the  atoms  from  the  origin  C^,  i.e. , . . . are  fixed; 

and  when  a ratio  P/n  has  been  decided  upon,  the  pair  of  positions  i which  must 
lie  on  the  circle  are  also  fixed,  together  with  the  directions  of  the  x and  z 

axes  and  the  angle  of  tilt  ^ = P/(nL).  Hence  the  co-ordinates  will  be  given  by  ex- 
pressions of  the  form 

2)  L.  Pauling,  R,  B.  Corey  and  H.  R.  Branson,  Proc.  Nat.  Acad.  Sci.  Wash. 

37,  205  (1951). 
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X-  = 0 

2^ 

= 0 

^1 

XQj  = CjOj  cos  (<j>  + 0^) 

= CjOj  sin  + 0j) 

*C|  “ SS  ®2^ 

= sin  + 0^) 

x^  = cos  ^ = L 

'c. 

= sin  ^ = p/n 

= projected  length  of  unit 


for  the  series  of  helices  with  an  upward  tilt,  while  in  the  other  series  the 
signs  of  the  angles  0^,  . . . will  be  reversed. 

All  these  co-ordinates  will  be  smooth  functions  of  and  hence  of  P/n, 
and  may  therefore  be  plotted  against  P/n  for  the  range  of  values  to  be  con- 
sidered. as  has  been  done  in  Fig.  2.  Values  of  the  co-ordinates  for  the  ratio 
p/n  = 1.48  lie  on  the  dotted  line  a 


Fig.  3.  Helix  projected  on  a 
plane  perpendicular 
to  the  helical  axis. 


= L + 1 + cos 
= L [p  + 0 + sin 
= (r  - 1)  P/n  . 


Second  stage.  Figure  3 shows  the 
helix  projected  on  a plane  perpendicular  to 
the  helical  axis.  The  vertical  plane  con- 
taining each  structural  unit  becomes  a line 
of  length  and  the  angle  between  two  such 
lines  is  the  helical  angle  ^ - 360/n,  and 
is  independent  of  the  pitch  P.  Since  these 
lines  are  equidistant  from  the  helical  axes, 
their  ends  C^,  C^,  . . . lie  on  the  circum- 

ference of  a circle  of  radius 

R = y L cosec  ^/2  . (1) 

The  co-ordinates  of  the  origin  of  the  rth  unit 
are  thus 

^ + cos  Z \jf  •¥  . . . 4 cos  (r  ' 2) 

\j/  + sin  Z \f/  -t  . . . + sin  {r  • Z)  J 


(2a) 

(Zb) 


•iipv  -4  - 


- 7 - 


The  co-ordinates  of  the  point  in  the  rth  unit  are  given  by 

% = *c 

*^r  r n 

Vp^  = YC^  + Xp^  sin  (r  - 1)  (3b) 

z = (r  - 1)  P/n  + z . (3c) 

‘'r  Pi 

Equations  (2a)  and  (2b)  can  be  rewritten  in  the  form 

= Lu^  . . . (4a) 

Vq  = ...  , (4b) 

r 

where  u^  and  are  the  unitary  co-ordinates  of  the  origins  of  successive  units, 
i.  e.  , 

u^  = fl  + cos  {jf  + cos  2\j/+  . . . cos  (r  - 2)  l/f]  (5a) 

v^  = fsin  [jf  + sin  Z\jf+  . . . sin  (r  - 2)  l/>]  (5b) 

These  functions  are  shown  in  Fig.  4 for  r = 2,  3,  4,  5,  6 over  the  range  n = 3.  0 


to  5.5. 


Equations  (3a)  and  (3b)  may  be  rewritten 


*P  ""  ^“r  ^ \ *p  • • • 

Yp  = Lv^  ®r  • • • ’ 

r Pi 

where  = cos  (r  - 1)  \jf  and  = sin  (r  - 1)  . Values  of  and  for  r = 2 

to  6 are  shown  in  Fig.  5 for  the  range  n = 3.  0 to  5.  5. 

Hence,  to  determine  atomic  co-ordinates  in  the  rth  unit  in  a helix  of  pitch 

P with  an  n-fold  screw:  (a)  the  co-ordinates  x and  z in  the  basic  unit  and  the 

Pi  Pi 

value  of  1^  for  the  requisite  value  of  P/n  are  obtained  from  Fig.  2;  (b)  the  values 
of  u^  and  v^  obtained  from  Fig.  4 for  the  chosen  value  of  n and  multiplied  by  L; 
(c)  the  values  of  A^  and  B^  are  obtained  from  Fig.  5 for  the  chosen  value  of  n 
and  multiplied  by  the  value  of  ;(d)  the  co-ordinates  are  written  down  from 
Eqs.  (6a),  (6b)  and  (3c). 


Since  the  co-ordinates  have  been  obtained  in  an  orthogonal  system,  bond 
lengths  and  bond  angles  can  be  determined  from  the  usual  expressions,  so  that 


Fig.  4.  Unitary  co-ordinates  of  origins  of  successive  structural  units  as  functions  of  n. 
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the  distance  ^ between  two  atoms  at  (x^y^z^)  and  (x^y^z^)  is  given  by 

®f.2  = + t^2  ‘ + t*2  ■ * 

and  the  angle  6 subtended  at  (x^y^z^)  by  the  atoms  at  (x^y^z^)  and  (x^y^Zj)  is 
given  by 


cos  6 = 


(*3  - *2^  ^*2  ■ *1^  ^^3  ■ ^2^  ^^2  ■ " *2^  ^®2  " *1^ 


°3.2‘  °2,1 


These  expressions  can  be  rewritten  in  various  ways  in  terms  of  the  values  of 
X,  y and  z given  by  Eqs.  (6a),  (6b)  and  (3c)  as  a specific  problem  may  require. 

For  some  purposes,  such  as  the  calculation  of  the  radial  distribution 
function,  it  may  be  preferable  to  express  the  atomic  positions  in  terms  of 
cylindrical  co-ordinates,  with  the  origin  of  the  co-ordinate  system  taken  on  the 
helical  axis. 

As  can  be  seen  in  Fig.  3,  all  the  atoms  of  a particular  type  lie  on  a circle 

of  a particular  radius,  and  it  can  be  seen 

from  Fig.  6 that  this  radius  p is  given  by 

p^  = R^+x^-  2RXp  cos  (90  - l/r/2) 

= + Xp  - 2RXp  sin  ( \^/2). 

But  R = ^ L cosec  ( '/'/2). 

p^  = ^ + x^  - Lx  . 

^ 4 sin  {[Jf/Z)  ^ P 

And  if  OC  is  taken  as  the  line  a = 0,  then 
a 

P 

The  foregoing  section  gives  analytical 
expressions  for  the  co-ordinates  of  the  atoms  in  the  helix.  It  can  be  seen  in 
Fig.  7,  however,  that  if,  in  addition  to  the  circle  of  diameter  drawn  in 

Fig.  1.  circles  having  diameters  ^1^1  ^1^1  drawn,  any 

choice  of  P/n,  which  fixes  the  points  (a)  and  (y)  results  in  a line 

4r  4T 

which  cuts  the  extra  circles  at  the  points  O^,  , etc. 


0 


Fig.  6.  Nomenclature  for 

polar  co-ordinates. 
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Hence  the  lengths  C^O*.  C^H*  etc.,  are  the  x co-ordinates,  and  OjO*, 
HjHj , etc.,  are  the  z co-ordinates  in  the  basic  unit.  Co-ordinates  of  atoms 
in  successive  units  can  be  obtained  either  by  utilizing  Figs.  4 and  5,  or 
graphically  by  constructing  the  projection  of  the  helix,  and  reading  the  co- 
ordinates directly  by  superposing  the  tracing  on  a sheet  of  graph  paper. 

However,  since  the  construction  of  such  a helix  involves  the  repetition 
several  times  of  selected  constants,  such  as  the  helical  angle,  it  is  particularly 
liable  to  cumulative  errors;  hence,  while  a graphical  construction  is  valuable  as 
an  aid  to  visualizing  the  general  feasibility  of  a postulated  helix,  it  does  not 
seem  suitable  for  obtaining  accurate  values  of  the  co-ordinates. 


